Structural and Magnetic Properties of Trigonal Iron 



S. Fox and H.J.F. Jansen 
Department of Physics, Oregon State University, Corvallis, OR 97331, USA 

(February 8, 2008) 

Abstract 

First principles calculations of the electronic structure of trigonal iron were 
performed using density function theory The results are used to predict lat- 
tice spacings, magnetic moments and elastic properties; these are in good 
agreement with experiment for both the bcc and fee structures. We find 
however, that in extracting these quantities great care must be taken in inter- 
preting numerical fits to the calculated total energies. In addition, the results 
for bulk iron give insight into the properties of thin iron films. Thin films 
grown on substrates with mismatched lattice constants often have non-cubic 
symmetry. If they are thicker than a few monolayers their electronic struc- 
ture is similar to a bulk material with an appropriately distorted geometry, 
as in our trigonal calculations. We recast our bulk results in terms of an iron 
film grown on the (111) surface of an fee substrate, and find the predicted 
strain energies and moments accurately reflect the trends for iron growth on 

a variety of substrates. 
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I. INTRODUCTION 



Thin transition metal films can be grown on a variety of substrates. The position of the 
atoms in the first film layer is largely determined by the surface geometry of the substrate. 
As more layers are added the lattice geometry of the film is determined by the forces at 
the substrate interface and by the inter-atomic forces within the film. As the film gets even 
thicker its geometry will eventually become that of the stable bulk material. Films with just 
a few layers give us the opportunity to study transition metals in geometries other than that 
of the bulk material, but with the additional complication of a interface with the substrate. 
Electronic structure calculations that include the effects of the substrate interface have 
revealed that other than the imposition of a non-bulk lattice geometry throughout the film, 
only the interface layer of the metal film is strongly affected by the substrate.S As a result 
bulk electronic structure calculations which ignore the details of the interface interactions 
but which use lattice constants characteristic of the substrate should give a good picture of 
the thin film electronic structure. 

Thin films of iron can be grown on substrates with a variety of lattice constants, so a 
range of calculations, corresponding to several different distortions of the bulk structure, 
are necessary to characterize these films. A distortion of a cubic unit cell along the (111) 
axis can be examined using a trigonal basis. The trigonal basis, characterized by three 
equal length vectors originating at the origin and the equal angle which separates each of 
them, contains the bcc, fee, and simple cubic structures as special cases. By varying the 
characteristic angle one can obtain these special geometries and the distortion along the 
(111) axis which connects them. Electronic structure calculations of distorted structures 
near the bcc and fee geometries will allow us to extrapolate elastic constants for these bulk 
geometries. More importantly in this context, these trigonal structures are exactly those 
which match the lattice provided by the (111) face of an fee substrate. So an electronic 
structure calculation of bulk iron with a trigonal geometry gives us insight into iron films 
grown on these substrates. Similarly, a tetragonal unit cell provides the correct distortion 
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to match the (100) face of a cubic substrate.B 

In this brief report we will give the results of electronic structure calculations for ferro- 
magnetic bulk iron over a range of trigonal geometries. The full-potential linearized aug- 
mented plane wave method (FLAPW) was used with Janak's parametrization of von Barth 
and Hedin's local density approximation to the exchange and correlation potential.i'S Some 
refinement of results from similar previous calculations for a tetragonal geometry will also 
be presented]! 



II. RESULTS FOR BULK MATERIALS 

Figure [I] shows the range of lattice geometries over which the energy calculations were 
performed. For this one atom per unit cell trigonal structure the geometry is parameterized 
in terms of the c/a ratio and the volume per unit cell. The natural log of the c/a ratio is 
used since it provides a more natural scaling. The volume is scaled by the experimental bcc 
volume (78.83 a.u.). The three special geometries, fee, simple cubic, and bcc, correspond to 
ln(c/a) values of 0.69, 0.0 and -0.69 respectively. Geometries at which electronic structure 
calculations were performed are indicated by a dot. Using the total energy result at each of 
these points we can numerically construct constant energy contours as shown in figure |l|. 

The most prominent features in this contour plot are the two local energy minimas cor- 
responding to an fee and a bcc structure. As has been observed in previous calculations the 
global minima is mispredicted; the fee structure lies about 4 mRyd low than the bcc struc- 
ture. At lower temperatures iron is actually a bcc crystal so we would expect a bcc structure 
to appear as the global minima for these calculations. This error is widely attributed to the 
inherent inaccuracy of the specific form of the local density approximation we have usedH 
One can travel between these two locally stable structures by varying the c/a ratio, which 
corresponds to distorting either cubic structure along the (111) axis. The energy barrier 
between bcc and fee structures for this distortion is approximately 70 mRyd/atom, signif- 
icantly larger than the 10 mRyd/atom found for the tetragonal distortion which connects 
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the bcc and fee structures.^ 

The data near the two minima can be used to obtain lattice constants, volumes and elastic 
properties of these two high symmetry structures. Ideally one would investigate regions 
extremely close to the energy minima where the energy varies quadratically with small 
distortions in the geometry. In this case a simple quadratic fit would perfectly model the 
behavior of the electronic energy around the minima allowing the various crystal properties 
to be extracted with great reliability. Unfortunately, because of the limited accuracy in the 
calculated energies, we need to look at structures further from the minima to ensure the 
energy differences among our structures are larger than the background numerical noise. As 
a result, we can no longer expect a simple quadratic fit to perfectly model variation in total 
energy among the structures near the minima. 

To account for this inherent inaccuracy the data were fit by a number of methods. Two 
dimensional fits of energy versus volume along the fee (ln(c/a) = 0.69) and bcc (ln(c/a) = 
—0.69) lines were performed with a simple quadratic model as well as with the Murnaghan 
and Birch-Murnaghan models.! Three-dimensional quadratic fits were made to points near 
the minima in the c/a versus V plane. The particular points to which each model was 
fitted were varied and the sensitivity of the final results to these variations was noted. 
In all cases 'good' (in terms of Chi-squared) fits were obtained. The variation in results 
between the different fits was used as an estimate of the uncertainty introduced by the use 
of imperfect fitting models. The lattice constants and volumes obtained were consistent for 
the various fitting techniques. The elastic constants showed some variation, especially for 
the fee minima. The elastic constants obtained from the three-dimensional quadratic fit were 
especially sensitive to the selection of fitted points. Tables | and |I| show some sample results 
from a variety of fits done near the bcc and fee minima. Considered carefully, the fitting 
errors are no greater than those known to be inherent in the local density approximation: 
2-3% in the lattice constants and 10-20% in the elastic constants. 

Combining these fit results and similar results from a refitting of previous calculations 
on tetragonally distorted iron, we obtained the predicted bulk ferromagnetic iron properties 
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shown in table |nT[ § Along with our results, table |T| gives experimental numbers for these 
properties along with the results from other density functional calculations.!'! The numbers 
for the bcc minima are in good agreement with the other theoretical predictions (also done 
under the LDA) consistently underestimating the lattice constants by a few percent, and 
over-estimating the elastic properties by a few tens of percent. For the fee data we present 
theoretical numbers from the so called 'low-spin' ferromagnetic state and experimental re- 
sults from gamma-iron, the high temperature fee phase. The discrepancy between theory 
and experiment is larger in this comparison, which is expected since the theoretical calcu- 
lations do not take temperature fully into account. In both the fee and bcc case the large 
values for the elastic constants can be largely attributed to the underestimate of the lattice 
constant. In general one finds in LDA calculations that the smaller the equilibrium volume, 
compared to the experimental value, the larger the bulk modulus. If we were to force the 
unit cell volume to be closer to the experimental value we could obtain elastic properties 
closer to those experimentally measured. 



III. RELATION TO THIN FILMS 

As was mentioned above, these bulk results can be related to thin iron films by recasting 
the data in terms of iron grown in perfect registry with the (111) surface of an fee substrate. 
In figure we show a contour plot of the same data as figure p], but now it has been 
parameterized in terms of a film geometry. The in-plane lattice spacing in the film will be 
fixed by the substrate since we assume the film is grown pseudomorphically. The interlayer 
spacing will take on whatever value minimizes the energy for the given in-plane spacing. If 
we know the lattice constant of our substrate we can immediately extrapolate from figure 
^| the interlayer spacing that minimizes the energy. Following this procedure for a range 
of substrate lattice constants we construct figure |3] which shows the energy of the film 
with optimal interlayer spacing as a function of the substrate lattice constant. The left 
and right minima in this figure correspond to unstrained fee, and unstrained bcc growth 
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respectively. The energy shift near these minima can be interpreted as the strain energy per 
layer associated with strained, pseudomorphic growth. 

As shown in figure |] we can, following a similar procedure, express the results of our 
bulk magnetic moment calculations as a function of substrate lattice constant, again with 
optimal interlayer spacing. In this figure we can see that unstrained fee films have small 
moments, but increasing the substrate lattice constant causes the moment to rise. In the 
region of bec growth the moment is not particularly sensitive to changes in the substrate. 
Keep in mind that these results ignore the interface effects which play a key role in the 
magnetic properties of such films. 

With this caveat in mind we can relate these figures to the observered properties of thin 
iron films for a variety of substrates. Copper, with an fee lattice constant of 6.82 a.u. should 
provide a reasonable match for pseudomorphic growth of iron films. Figures 0, |3|, and |] 
predict iron grown on copper (111) is a strained fee structure with an interlayer spacing of 
around 3.7±0.2 a.u. and a moment of 1.0±0.25. Recent studies showed strained fee growth 
up to five layers with interlayer spacings of 3.9 a.u. and moments of lfj,B/&tom consistent 
with our result Less well matched substrates have been studied as well. Using Pd , 
which has a lattice constant of 7.36 substrate should create a strain energy per 

layer nearly twice that of copper according to figure |[ Not surprisingly it was found that 
less than two layers of iron could be grown pseudomorphically on Pd before islands of bec 
iron started to form.0 Substrates which would provide even larger strains than Pd, namely 
Ag and Al do not appear to give clean pseudomorphic growth as expected from figure ^.0 

IV. CONCLUSIONS 

In this paper we have examined a range of trigonal iron structures using density functional 
theory. We find energy minima corresponding to the expected bec and fee geometries. 
While carefully considering the inherent inaccuracy in the fitting process we can extract 
information about the properties of these bulk states which are in good agreement with 
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experimental results. We also recast the bulk results to reflect the geometries associated 
with the pseudomorphic growth of thin iron films. Although this approach ignores the 
interface and surface effects it still gives meaningful insight into the growth of iron films on 
a range of substrates. The predicted film properties are consistent with the growth observed 
for a variety of substrates. 
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FIGURES 

FIG. 1. Contours of constant energy (mRyd/atom) of ferromagnetic trigonal iron as a function 
of volume and c/a ratio. Dots indicate geometries at which electronic structure calculations were 
performed. All energies are offset by 2541 mRyd. 

FIG. 2. Contours of constant energy (mRyd/atom) of a ferromagnetic iron film as a function 
of fee substrate lattice constant and film interlayer spacing. Dots indicate geometries at which 
electronic structure calculations were performed. Contours far from data points are unreliable. All 
energies are offset by 2541 mRyd. 

FIG. 3. Energy (mRyd/atom) of a thin iron film with optimal interlayer spacing as a function 
of the fee substrate lattice constant. All energies are offset by 2541 mRyd. A line is provided to 
guide the eye. 

FIG. 4. Magnetic moment (fiB /atom) of a thin iron film with optimal interlayer spacing as a 
function of the fee substrate lattice constant. 
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TABLES 





TABLE I. Properties of bcc iron as predicted by 




c/a 


Volume (a.u.) 


Energy (mRyd.) 


Bulk Modulus (Mbar 


) r it ting Met no a 


U.O 


73.4 


-0.103 


2.25 


parabola to 4 points 


U.O 


71.9 


-0.102 


1.77 


Murnaghan near minima 


U.O 


73.3 


-0.103 


2.30 


Birch-Murnaghan near minima 


0.5 


71.1 


-0.103 


2.19 


Murnaghan over full range of V 


0.5 


71.2 


-0.104 


2.33 


Birch-Murnaghan over full range of V 


0.5084 


73.2 


-0.101 


2.11 


3-d quadratic with minimum error 


0.5092 


73.1 


-0.101 


2.14 


3-d quadratic with average error 


0.5101 


73.1 


-0.102 


2.22 


3-d quadratic with average error 



TABLE II. Properties of fee iron as predicted by a variety of fitting models. 



c/a Volume (a.u.) Energy (mRyd.) Bulk Modulus (Mbar) Fitting Method 



2.0 65.8 


-0.105 


3.19 


parabola to 4 points 


2.0 65.9 


-0.109 


3.09 


Murnaghan near minima 


2.0 65.7 


-0.105 


3.16 


Birch-Murnaghan near minima 


2.0 66.0 


-0.103 


2.20 


Murnaghan over full range of V 


2.0 65.6 


-0.104 


2.57 


Birch-Murnaghan over full range of V 


2.023 67.2 


-0.107 


3.66 


3-d quadratic with minimum error 


2.006 66.9 


-0.111 


3.74 


3-d quadratic with average error 


2.014 66.9 


-0.103 


2.30 


3-d quadratic with average error 
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TABLE III. Properties of bcc and fee iron. Results from this study and other density functional 
calculations, as well as experimental results for a-iron and 7-iron. All elastic properties are in Mbar. 



lattice constant (a.u.) moment (/uB) Bulk Modulus c\\ cyi C44 

bcc 

theory 5.24 ±0.1 2.0 ± 0.05 2.2 ±0.15 2.9 ±1.1 1.8 ±0.6 1.5 ±0.3 

theory a 5.21 2.08 2.66 

exper. b 5.406 2.12 1.68 2.33 1.355 1.178 

fee 

theory 6.43 ± 0.1 0.0 2.6 ± 0.4 5.4 ±1.0 1.2 ±0.3 2.6 ± 0.6 

theory c 6.37 0.0 3.9 

exper. d 6.947 0.0 1.32 1.54 1.22 0.77 

a from Ref. | 
b a-iron from Ref. [?] 
c Low-Spin results from Ref. || 
d 7-iron from Ref. 
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Trigonal Iron — Constant Energy Contours 
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